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Abstract 



The IceCube Observatory is a kilometer-cube neutrino telescope under construction 
at the South Pole and planned to be completed in early 2011. When completed it will 
consist of 5,160 Digital Optical Modules (DOMs) which detect Cherenkov radiation 
from the charged particles produced in neutrino interactions and by cosmic ray initiated 
atmospheric showers. IceCube construction is currently 90% complete. A selection of 
the most recent scientific results are shown here. The measurement of the anisotropy 



C3 , in arrival direction of galactic cosmic rays will also be presented and discussed. 



1 Introduction 



The IceCube Observatory is a km 3 neutrino telescope designed to detect astrophysical 
neutrinos with energy above 100 GeV. IceCube observes the Cherenkov radiation from 
charged particles produced in neutrino interactions. 

The quest for understanding the mechanisms that shape the high energy Universe is 
taking many paths. Gamma ray astronomy is providing a series of prolific experimental 
observations, such the detection of TeV 7 rays from point-like and extended sources, along 
with their correlation to observations at other wavelengths. These observations hold the 
clues about the origin of cosmic rays and the possible connection to shock acceleration in 
Supernova Remnants (SNR), Active Galactic Nuclei (AGN) or Gamma Ray Bursts (GRB). 
Supernovse are believed to be the sources of galactic cosmic rays, nevertheless the 7 ray 
observations from SNR still do not provide us with a definite and direct evidence of proton 
acceleration. The competing inverse Compton scattering of directly accelerated electrons 
may significantly contribute to the observed 7 ray fluxes, provided that the magnetic field 
in the acceleration region does not exceed 10 /iG pQ. 

Ultra High Energy Cosmic Rays (UHECR) astronomy, has the potential to hold the 
key to a breakthrough in astroparticle physics. The identification of sources of cosmic rays 



could provide a unique opportunity to probe the hadronic acceleration models currently 
hypothesized. On the other hand, cosmic ray astronomy is only possible at energies in 
excess of 10 19 eV, where the cosmic rays are believed to be extragalactic and point back to 
their sources. TeV 7 rays from those sources are likely absorbed during their propagation 
between the source and the observer : at ~10 TeV 7 rays have a propagation length 
of about 100 Mpc, while at ~100 GeV 7 rays can propagate much deeper through the 
Universe. 

If the extra-galactic sources of UHECR are the same as the sources of 7 rays, then 
hadronic acceleration is the underlying mechanism and high energy neutrinos are produced 
by charged pion decays as well. Neutrinos would provide an unambiguous evidence for 
hadronic acceleration in both galactic and extragalactic sources, and are the ideal cosmic 
messengers, since they can propagate through the Universe undeflected and with practically 
no absorption. But the same reason that makes neutrinos ideal messengers makes them 
also difficult to detect. 

The discovery of the anisotropy in arrival direction of the galactic cosmic rays has 
also triggered particular attention recently. The origin of galactic cosmic ray anisotropy is 
still unknown. The structure of the local interstellar magnetic field within 1 pc is likely 
to have an important role in shaping the large angular scale features of the observed 
anisotropy. Nevertheless it is possible to argue that the anisotropy might be originated 
by a combination of astrophysical phenomena, such as the distribution of nearby recent 
supernova explosions [2j. The observation of galactic cosmic ray anisotropy at different 
energy and angular scales has, therefore, the potential to reveal the connection between 
cosmic rays and shock acceleration in supernovae. 

At the same time, there seems to be clear observational evidence for the existence of dark 
matter in the Universe, even if its nature remains unknown. A variety of models predict the 
existence of a class of non-relativistic particles called Weakly Interacting Massive Particles 
(WIMP). These particles could be gravitationally condensed within dense regions of matter 
(such as the Sun or the galactic halo) and could provide a visible source for indirect 
detection via neutrino generation through annihilation. Neutrino telescopes are powerful 
tools to indirectly test spin-dependent WIMP-nucleon scattering cross section, provided 
the models for matter distribution and WIMP annihilation rate are taken into account. 

In section $2] the IceCube Observatory apparatus functionality and calibration are de- 
scribed. Selected physics analyses results are summarized in $3] : the determination of the 
atmospheric muon neutrino energy spectrum (£ )3.ip . the search for astrophysical neutrinos 
from diffused and point sources, and from Gamma Ray Bursts ( §3.2p . the indirect search 
for dark matter (£ |3.3p . and the anisotropy in cosmic rays arrival direction ( §3.4p . 



2 The IceCube Observatory 

The IceCube Observatory (see Fig. [TJ currently consists of 4,740 DOMs deployed in 79 
vertical strings (60 DOMs per string) between 1,450 m and 2,450 m depth below the 
Geographic South Pole. At the beginning of 2011 IceCube will be completed with 86 
strings and 5,160 DOMs. The surface array IceTop with 81 stations, each consisting of 
two tanks with frozen clean water with two DOMs each, will provide the measurement of 
the spectrum and mass composition of cosmic rays at the knee and up to about 10 18 eV. 
The Deep Core sub-array, consisting of 6 densely instrumented strings and located at the 
bottom-center of IceCube, is capable of pushing the neutrino energy threshold to about 10 
GeV. The surrounding IceCube instrumented volume can be used to veto the background 
of cosmic ray induced through-going muon bundles to enhance the detection of down-going 
neutrinos within the Deep Core volume. The veto rejection power can reach 10 5 . 
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Figure 1: The IceCube Observatory. Currently it consists of 79 strings and 4,740 DOMs. 
In 2011 it will be completed with 86 strings and 5,160 DOMs. The shaded region near 
the center of the array is the Deep Core dense sub-array, and the one on the right is the 
AMANDA neutrino telescope, de-comissioned in May 2009. 



AMANDA was the first and largest neutrino telescope before the construction of Ice- 
Cube. Using analog technology it significantly contributed to the advance of neutrino 



astrophysics searches and it was de-commissioned in May 2009. 

The basic detection component of IceCube is the DOM : it hosts a 10-inch Hamamatsu 
photomultiplier tube (PMT) and its own data acquisition circuitry enclosed in a pressure- 
resistant glass sphere, making it an autonomous data collection unit. The DOMs detect, 
digitize and timestamp the signals from optical Cherenkov radiation photons. Their main 
board data acquisition (DAQ) is connected to the central DAQ in the IceCube Laboratory 
at the surface, where the global trigger is determined [6]. 

The detector calibration is one of the major efforts aimed at characterizing its response 
and to reduce systematic uncertainties at the physics analysis level. Each PMT is tested 
in order to characterize its response and to measure the voltage yielding a specific gain [7j. 
In the operating neutrino telescope the current gain is about 10 7 and the corresponding 
dark noise rate is about 500 Hz. Time calibration is maintained throughout the array by 
regular transmission to the DOMs of precisely timed analog signals, synchronized to a cen- 
tral GPS-disciplined clock. This procedure has a resolution of less than 2 nsec. The LEDs 
on the flasher boards instrumented in the DOMs, are used to measure the photo-electron 
(p.e.) transit time in the PMT for the reception of large light pulses between neighboring 
DOMs. This delay time is given by light travel time from the emitter to the receiver, by 
light scattering in ice and by the electronics signal processing. The RMS of this delay is 
also less than 2 nsec. Waveform sampling amplitude and time binning calibration is peri- 
odically performed in each DOM and used to extract the number of detected p.e. with an 
uncertainty of less than 10%. Higher level calibrations are meant to correlate the number 
of detected p.e. to the energy of physics events that trigger the detector. Instrumented 
devices, such as the flasher boards, are used to illuminate the detector with 400 nm wave- 
length photons (corresponding to the wavelength yielding the highest detection sensitivity), 
simulating a real electron-neutrino interaction, or cascade, inside the detector. A complete 
Monte Carlo simulation chain is used to relate the known number of injected photons with 
the energy scale of the artificial cascade. The energy resolution depends on the event topol- 
ogy (track-like versus cascade-like) and on its containment inside the instrumented volume. 
Monte Carlo simulations provide the necessary fluctuations implicit by the topology and 
containment of the physics events. The ice optical properties are the most fundamental 
calibration determination that allows us to know how photons propagate through the ice 
and, therefore, how to relate the number of detected p.e. to the energy of the physics 
events. Due to the antarctic glaciological history, the optical properties depend on depth 
and they have been measured in the past using AMANDA in-situ calibration lasers [8], in 
the depth range between 1,400 m and 2,000 m. The optical properties down to 2,450 m, 
the depth of IceCube instrumentation, are extrapolated from ice core observations in other 
location of the antarctic continent, and a new campaign of extended in-situ measurements 
is currently being carried out. 



3 Physics Results 

If the DOMs that detect Cherenkov photons satisfy specific trigger conditions, an event is 
formed and recorded by the surface DAQ. An on-line data filtering at the South Pole reduces 
the event volume to about 10% of the trigger rate, based on a series of filter algorithms 
aimed to select events based on directionality, topology and energy. The filter allows us to 
transfer data via satellite to the northern hemisphere for prompt physics analyses. 

3.1 Atmospheric neutrinos 

99.999% of the events that trigger IceCube are muons produced by the impact of primary 
cosmic rays in the atmosphere. Only a small fraction of the detected events (~10~ 5 ) are 
muon events produced by atmospheric neutrinos. In order to reject all the down-ward 
muon bundle background only up-ward events are generally selected, assuming specific 
event selections provide well reconstructed events. In the 40 string instrumented IceCube 
(IceCube-40) about 30-40% of the up-ward events survive the selection, with a background 
contamination of less than 1% |9|. 
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Figure 2: The preliminary unfolded spectrum of atmospheric muon neutrinos {y^+Va) from 
the IceCube-40 detector [9], compared with previously measurements : the Frejus result 
|10j . upper and lower bands from Super-Kamiokande |llj . the forward-folding result from 
AMANDA p2] and the unfolded spectrum from AMANDA [IS]. The AMANDA unfolding 
analysis was a measurement of the zenith-averaged flux from 100° to 180° and the present 
analysis (IC40 unfolding) is a measurement of the zenith-averaged flux from 97° to 180°. 



The energy estimation resolution of these atmospheric v„ induced events is of the order 
of 0.3 in Log of neutrino energy, and a regularized unfolding technique is used to determine 



the energy spectrum. Fig. [2] shows the preliminary unfolded energy spectrum of the 
17,682 atmospheric neutrinos detected by IceCube-40 between zenith angle of 97° and 
180°. The figure also shows other measurements performed by other experiments, including 
AMANDA |12|I13]. IceCube has detected the highest energy atmospheric neutrinos (about 
250 TeV where a significant fraction of neutrinos are expected to arise from the decay of 
heavy mesons with charm quarks). IceCube allowed us to extend the global measured 
spectrum up to 6 orders of magnitudes in energy. For the first time the precision of this 
measurement is providing a powerful tool to test the high energy hadronic interaction 
models that govern our present knowledge of the cosmic ray induced extensive air showers. 

3.2 Search for astrophysical neutrinos 

Atmospheric neutrinos represent an irreducible background for the search of high energy 
astrophysical neutrinos. If hadronic acceleration is the underlying source of high energy 
cosmic rays and 7 rays observations, we expect that unresolved sources of cosmic rays over 
cosmological times have also produced enough neutrinos to be detected as a diffuse flux. 
Since shock acceleration is expected to provide an ~E~ 2 energy spectrum, harder than 
the ~E~ 3 - 7 of the atmospheric neutrinos, the diffuse flux is expected to dominate at high 
energy. 
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Figure 3: Current status on the search of diffuse astrophysical v^ + v^. The atmospheric 
neutrino background as measured by IceCube-40 and AMANDA are shown, along with the 
predictions to the highest possible energies [14l fl~5| [IBj I17j . The limit on an E~ 2 flux by 
AMANDA [TS] and the preliminary sensitivity by IceCube-40 are shown. We predict that 
IceCube-40 sensitivity is below the Waxman & Bahcall diffuse neutrino bound [19]. The 
graph shows other models of astrophysical muon neutrinos as well [20\ |2"T| |2"21 12"3"| 124] (An 
E~ 2 spectrum is a horizontal line in this graph). 



Fig. [3] shows the AMANDA experimental upper limit and the preliminary IceCube-40 



sensitivity for an E~ 2 diffuse muon neutrino spectrum (y^ + tv). One year of IceCube-40 
is about 5 times more sensitive than 3 years of AMANDA and its sensitivity is lower than 
the Waxman-Bahcall neutrino bound. This means that IceCube is potentially approaching 
the discovery of the origin of cosmic rays. 

In the Ultra High Energy range (i.e. above ~10 6 GeV or UHE) IceCube is placing 
upper limits that are still more than about one order of magnitude above the predicted 
flux of neutrinos from UHECR interactions on the microwave photons (or GZK neutrinos 
|25|). The complete IceCube Observatory might be able to reach the discovery level within 
the next 5-8 years. 

If the observed 7 rays from galactic and extra-galactic point and extended sources are 
from neutral pion decays in hadronic acceleration sites or from cosmic ray interaction with 
molecular clouds, the charged pions could produce enough neutrinos that can be detected 
by a km 3 neutrino telescope. 
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Figure 4: On the left: sensitivity (90% CL) for a full-sky search of steady point sources 
of muon neutrinos with an E -2 energy spectrum as a function of declination angle (5 < 
0° for down-ward events, or the southern hemisphere, and 5 >0° for up-ward events, 
or the northern hemisphere). The IceCube-22 sensitivity for the full-sky search |26|. I27j 
is about the same as the IceCube-40 discovery potential (5a discovery for 50% of the 
trials), and about a factor of three larger than the IceCube-40 sensitivity (90% CL). The 
squares show the IceCube-40 90% upper limits for selected sources, and the black line 
represents the projected IceCube-86 sensitivity. On the right: preliminary sky map in 
equatorial coordinates of statistical significance (p-value) from the search for steady point 
sources of high energy muon neutrinos in IceCube-40. The observation is extended to 
the southern hemisphere by reducing muon background five orders of magnitudes with a 
zenith-dependent energy selection (>100's TeV). 



Fig. 0] shows, on the left, the sensitivity (90% CL) of IceCube for the full-sky search 
of steady point sources of E~ 2 muon neutrinos as a function of declination. The exten- 
sion of the point source search to the southern hemisphere is made possible by rejecting 
background events by five orders of magnitude with a high energy event selection. This 
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makes the southern hemisphere still dominated by high energy muon bundles and it yields 
a poorer neutrino detection sensitivity because of the high energy selection. Neverthe- 
less this opens IceCube to a full-sky coverage and provides a coverage complement to the 
neutrino telescopes in the Mediterranean. On the right of Fig. [5] is the sky-map of sta- 
tistical significance from the full sky search of IceCube-40. No significant localized excess 
is observed. The sensitivity is to be interpreted as the median upper limit we expect to 
observe from individual sources across the sky. If we test specific sources (see left panel of 
Fig. [5]) we see that the full IceCube (about twice as sensitive as IceCube-40) will be able 
to discover neutrinos from individual point sources in about 3-5 years, depending on the 
location in the sky. 
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Figure 5: On the left: sensitivity (90% CL) and discovery potential (5<r discovery for 50% 
of the trials) in IceCube-40 for the Crab [22], MGROJ1852+01 [25] and 3C279 [30 J. On 
the right: preliminary upper limits (90% CL) for the searches of neutrinos in coincidence 
with Gamma Ray Bursts in AMANDA [21], IceCube-22 [32] and IceCube-40. 



The search for neutrinos from transient galactic and extra-galactic sources is also being 
pursued. In particular, the right panel in Fig. [5] shows the upper limits (90% CL) for the 
model-dependent search of prompt neutrinos from GRBs in the northern hemisphere with 
AMANDA [21], IceCube-22 [32] and IceCube-40. For each detector configuration, the list 
of GRBs detected during the corresponding physics runs were collected and the predicted 
neutrino flux calculated based on the 7 ray spectrum [33J. The corresponding average 
neutrino spectrum was used to search for neutrinos detected within the so-called T$q time 
window (i.e. the time in which 5% to 95% of the fluence is recorded). The right panel of 
Fig. [5] also shows the Waxman-Bahcall (WB) predicted average spectrum from GRBs [34j 
and the average GRB spectrum corresponding to the 2009-2008 time period of IceCube-40 
physics runs. The preliminary IceCube-40 upper limit is below the WB spectrum, which 
indicates that IceCube is becoming very sensitive and could potentially discover neutrinos 
in coincidence with GRBs within the next few years. 



3.3 Search for dark matter 

Non-baryonic cold dark matter in the form of weakly interacting massive particles ( WIMPs) 
is one of the most promising solutions to the dark matter problem [35J . The minimal super- 
symmetric extension of the Standard Model (MSSM) provides a natural WIMP candidate 
in the lightest neutralino Xi [36J . This particle is weakly interacting only and, assuming 
R-parity conservation, is stable and can therefore survive today as a relic from the Big 
Bang. A wide range of neutralino masses, rm-.o, from 46 GeV [37] to a few tens of TeV 
|38j is compatible with observations and accelerator-based measurements. Within these 
bounds it is possible to construct models where the neutralino provides the needed relic 
dark matter density. 
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Figure 6: On the left: upper limits (90% CL) on the muon flux from neutralino annihila- 
tions in the Sun for the soft (66) and hard (W + W~) annihilation channels, adjusted for 
systematic effects, as a function of neutralino mass. A muon energy threshold of lGeV was 
used when calculating the flux. Also shown are the limits from MACRO [31], Super-K |42j . 
AMANDA [33], and IceCube-22 [44] merged to AMANDA at the low energy end. On the 
right: Upper limits (90% CL) on the spin dependent neutralino-proton cross-section a 
for the soft (66) and hard (W + W~) annihilation channels, adjusted for systematic effects, 
as a function of neutralino mass. Also shown are the limits from CDMS |39j, COUPP |45j, 
KIMS [3B] and Super-K [32]. The shaded area represents MSSM models not disfavored by 
direct searches [39j[30] based on a SI . 



Relic neutralinos in the galactic halo may be gravitationally attracted by the Sun 
and accumulate in its center, where they can annihilate each other and produce standard 
model particles, such as neutrinos. This provides an indirect channel detection of this type 
of dark matter, provided the WIMP density and velocity distribution and the neutralino 
annihilation rate models are taken into account. The left panel of Fig. [6] shows the upper 
limits (90% CL) on the muon flux for IceCube-22 [33] merged to the AMANDA upper 
limit at the low energy end [43J, along with other indirect observations. The limits on 
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the annihilation rate can be converted into limits on the spin-dependent a and spin- 
independent a neutralino-proton cross-section (as shown on the right panel of Fig. [U]). 
This conversion allows us a comparison with the direct search experiments. Since capture 
in the Sun is dominated by a SD , indirect searches are expected to be competitive in setting 
limits on this quantity. Assuming equilibrium between the capture and annihilation rates 
in the Sun, the annihilation rate is directly proportional to the cross-section. Fig. [6] 
also shows the predicted sensitivity (90% CL) of the combined IceCube-86 and the Deep 
Core dense instrumentation that allows us to significantly lower the energy threshold, and 
consequently increase sensitivity to low neutralino mass. In indirect searches WIMPs would 
accumulate in the Sun over a long period and therefore sample dark matter densities in a 
large volume of the galactic halo. This progressive gravitational accumulation is sensitive 
to low WIMP velocities, while direct detection recoil experiments are more sensitive to 
higher velocities, making indirect searches a good complement to the direct ones. 

Neutrino telescopes can also test the dark matter self-annihilation cross section < gav > 
(averaged over the dark matter velocity distribution), making them complementary to 7 
ray measurements. If the lepton excess observed by Fermi |47j . H.E.S.S. [H] and PAMELA 
[49] is interpreted as a dark matter self-annihilation signal in the galactic dark matter halo 
|50| . the leptophilic dark matter in the TeV mass range provides the most compatible 
model. Since the dark matter halo column density is larger toward the direction of the 
Galactic Center (GC), neutrinos from WIMP self-annihilation in the halo are expected to 
have a large angular scale anisotropy with an excess in the direction of the GC (see left 
panel of Fig. [7|). The dark matter density distribution in the Milky Way has different 
shapes depending on the model [51]. The expected neutrino flux from dark matter self- 
annihilation is proportional to the square of the dark matter density integrated along the 
line of sight for a given angular distance from the Galactic center. The differential neutrino 
flux for a WIMP of mass m x depends on the halo density profile, the neutrino production 
multiplicity, and on self-annihilation cross section < gav > [52]. The search for an excess 
of neutrinos in the direction of the GC for different annihilation channels (assuming 100% 
branching ratio for each of them), allows us to probe the allowed range of < ga v > for the 
corresponding channels in a model- independent manner (see right panel of Fig. [7]), and 
provide a direct comparison to similar results from 7 ray experiments. 

IceCube's reach can be significantly improved by directly looking at the GC (visible 
from the southern hemisphere). This search was performed with the IceCube 40-string 
dataset, by using down- ward neutrinos interacting inside the detector volume, and the 
preliminary upper limits (90% CL) on the self-annihilation cross-section are shown on the 
right of Fig [7] While such an analysis is already able to put significantly better constraints, 
a large scale anisotropy would provide a more distinct discovery signal. 
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Figure 7: On the left: relative expected neutrino flux from dark matter self-annihilations in 
the Milky Way halo on the northern celestial hemisphere (in equatorial coordinate). The 
largest flux is expected at right ascension a closest to the GC (Aa = 0). Dashed lines 
indicate circles around the Galactic center, while the solid lines show the definition of on 
and off-source region on the northern hemisphere. The on-source region is centered around 
Aa = 0, while the off-source region is rotated by 180° in a. On the right: preliminary 
upper limits (90% CL) on the self-annihilation cross-section < a^v > as function of the 
WIMP mass m x , for the search of an excess from the GC in the northern hemisphere with 
IceCube-22 |51j and the search for down-ward neutrinos from the GC interacting inside 
the IceCube-40 array, for different annihilation channels. The bands on the IceCube-22 
sensitivity curves account for different halo density profiles. The dashed line at the bottom 
labeled natural scale is for dark matter candidates consistent with being a thermal relic. 
The black dashed line is the unitarity bound. 



3.4 Cosmic ray anisotropy 

Galactic cosmic rays are found to have an energy dependent large angular scale anisotropy 
in arrival direction distribution with amplitude of about 10 -4 — 10~ 3 . The first compre- 
hensive observation of such an anisotropy was provided by a network of muon telescopes 
sensitive to sub-TeV cosmic ray energies and located at different latitudes [531 El] ■ More 
recently, an anisotropy was also observed in the multi-TeV range by the Tibet AS7 array 
[55], ARGO YBJ [56], Super-K [57] , and by MILAGRO [58]. And the first observation in 
the southern hemisphere is being reported by IceCube [59] for a median cosmic ray energy 
of about 20 TeV. 

The left panel of Fig. [8] shows the relative intensity in arrival direction of the cosmic 
rays, obtained by normalizing each ~3° declination band independently. On the top is 
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Figure 8: On the left: sky-map of the relative intensity in arrival direction of cosmic rays for 
IceCube-22 (on the top) [59], and preliminary sky-map of the relative intensity for IceCube- 
40 (on the bottom), in equatorial coordinates. For a better visual effect, a 3° smoothing has 
been applied to the maps. On the right: preliminary modulation of the relative intensity in 
arrival direction of cosmic rays for IceCube-40, projected in right ascension (black symbols), 
in the right ascension with respect to the Sun's location (green symbols), and in pseudo- 
right ascension, i.e. corresponding to anti-sidereal time (red symbols). 



the relative intensity map obtained from the 4.6 billion events collected by IceCube-22 [59J 
and on the bottom the preliminary map obtained from the 12 billion events collected by 
IceCube-40. The two maps show the same anisotropy features and they both appear to be 
a continuation of the observed modulation in the northern hemisphere. The right panel 
of Fig. [8] shows the preliminary modulation of the relative intensity in arrival direction of 
the cosmic rays projected into right ascension for IceCube-40 (black symbols). In order 
to verify whether the observed sidereal anisotropy (i.e. in equatorial coordinates), has, 
in one full year, some spurious modulation derived from the interference between possible 
yearly-modulated daily variations, the same analysis was performed using the anti-sidereal 
time frame (a non-physical time defined by switching the sign of the transformation from 
universal to sidereal time) [60]. The real feature in the sidereal time is expected to be 
scrambled in the anti-sidereal time. The anti-sidereal modulation (shown on the right of 
Fig. [8] in red symbols) appears to be relatively flat with an amplitude that is the same order 
of magnitude of statistical errors, suggesting that no significant spurious effect is present. If 
the relative intensity is measured, over one full year, as a function of the angular distance 
from the Sun in right ascension, we expect an excess in the direction of motion of the 
Earth around the Sun (at ~ 270°) and a minimum in the opposite direction. This is what 
is observed (see the green symbols on the right of Fig. [8]). 

With the same techniques used in 7 ray detection, it is possible to estimate the event 
intensity sky-map (the background) with a pre-defined angular scale averaging, and deter- 
mine the residual by subtracting the background from the actual map. With this method 
MILAGRO, in an attempt to estimate the background without 7-hadron separation, dis- 
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Figure 9: Preliminary statistical significance sky- map of the medium scale anisotropy in 
arrival direction of the cosmic rays for IceCube-40, combined with the MILAGRO northern 
hemisphere sky-map [61J. In this map only anisotropy scales that are smaller than ~30° are 
visible. Note the different color scales for the statistical significance in the two hemispheres. 



covered two significant localized regions of cosmic rays |61j . also observed by ARGO YBJ 
|62j . The same medium-scale anisotropy measurement was performed with IceCube for 
the first time in the southern hemisphere and the combined MILAGRO-IceCube-40 signif- 
icance sky-map is shown in Fig. [9l where only the anisotropy features with angular scale 
smaller than ~30° are visible. The different event statistics between MILAGRO (with 220 
billion events with median energy of ~ 1 TeV) and IceCube-40 (with 12 billion events with 
median energy of ~ 20 TeV) does not allow the comparison of the two hemispheres on 
a statistical base. Nevertheless, there seems to be some indication that the small scale 
features observed in the two hemispheres might be part of a larger scale structure. 

The origin of the galactic cosmic ray anisotropy is still unknown. However there might 
be multiple superimposed causes, depending on the cosmic ray energy and the angular 
scale of the anisotropy. The large scale structure in the 10-100 TeV range might be a 
local fluctuation caused by a nearby supernova (within 1,000 pc) that exploded within 
the last 100,000 years or so [2]. On the other hand, the structure of the local Inter-Stellar 
Medium (ISM) magnetic field well within 1 pc might likely have an important role [63]. The 
strongest and most localized of the MILAGRO excess regions has triggered astrophysical 
interpretations, invoking Geminga pulsar as a possible source |64}I65| . A strong anisotropy 
of the Magneto-Hydro Dynamic turbulence in the ISM, could cause a superposition of the 
large scale anisotropy (perhaps generated by a nearby SNR) with a beam of cosmic rays 
focused along the local magnetic field direction, depending on the turbulence scale [66J. 
However, the localized nature of the hottest MILAGRO excess region suggests a local origin. 
Neutron monitor data also seem to indicate that the broad excess toward the direction of 



13 



the heliotaiU (the so-called tail- in excess), which includes the MILAGRO localized excess 
regions, is likely generated within the heliotail itself [67]. In particular the tail-in excess and 
its small angular scale structure, is suggestive of acceleration via magnetic reconnection 
in the solar magnetotail. Reconnection is generated by magnetic polarity reversals due to 
the 11-year solar cycles compressed by the solar wind in the magnetotail. The maximum 
energy of protons that can be accelerated through this process is estimated to be about 10 
TeV. Up to this energy scale a localized excess might be observable in the direction of the 
acceleration sites [68] . This is the energy at which MILAGRO observes a cut-off for the 
localized regions. 
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